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Abstract. Recent optical measurements of deeply underdoped cuprates have revealed that a coherent 
Drude response persists well below the end of the superconducting dome. In addition, no large increase 
in optical effective mass has been observed, even at dopings as low as 1%. We show that this behavior is 
consistent with the resonating valence bond spin-liquid model proposed by Yang, Rice, and Zhang. In this 
model, the overall reduction in optical conductivity in the approach to the Mott insulating state is caused 
not by an increase in effective mass, but by a Gutzwiller factor, which describes decreased coherence due 
to correlations, and by a shrinking of the Fermi surface, which decreases the number of available charge 
carriers. We also show that in this model, the pseudogap does not modify the low-temperature, low- 
frequency behavior, though the magnitude of the conductivity is greatly reduced by the Gutzwiller factor. 
Similarly, the profile of the temperature dependence of the microwave conductivity is largely unchanged in 
shape, but the Gutzwiller factor is essential in understanding the observed difference in magnitude between 
ortho-I and -II YBa2Cu3 0j;. 

PACS. 74.72.-h cuprate superconductors - 74.25.Gz optical properties - 74.72.Kf pseudogap regime 



1 Introduction 

As formulated by Yang, Rice, and Zhang (YRZ) [1], the 
resonating valence bond model [2] of the superconduct- 
ing underdoped cuprates has been remarkably succcssbil 
in providing an understanding of properties previously 
considered anomalous [3-9]. An essential element of the 
model is the formation of a pseudogap that grows in mag- 
nitude with decreased doping and reconstructs the Fermi 
surface into Luttinger hole and electron pockets, which 
replace the usual large Fermi surface of Fermi liquid the- 
ory. This extension of BCS theory is quite diflFerent from 
the many others that were required to explain the prop- 
erties of conventional and even optimally and overdoped 
cuprates. Those extensions include anisotropy [10], energy 
dependence in the electronic density of states [11], inelastic 
scattering [12,13], strong coupling effects [14], and d-wave 
gap symmetry with possibly higher harmonics [15,16]. 

In this paper, we use the YRZ model to explain anoma- 
lous optical conductivity data. W.J. Padilla et al. [17] 
and Y.S. Lee et al. [18] recently measured the optical re- 
sponse in the deeply underdoped region of the high-Tg 
cuprate phase diagram. In both La2-xSra,Cu04 (LSCO) 
and YBa2Cu30j, (YBCO), they found that a coherent 
Drude response, characteristic of metallicity, persisted well 
below the end of the superconducting dome in the ap- 
proach to the antiferromagnetic Mott instulating state, 



even for doping values as low as x = 1% in LSCO and 
y = 6.3 in YBCO. They also noted that, remarkably, the 

optical effective mass of the charge carriers remains nearly 
constant, reaching at most a few times the bare electron 
mass. Hence, they concluded that the; approach to the in- 
sulating state must be due to a reduction of the number 
of charge carriers rather than a decrease in carrier mobil- 
ity. They suggested exotic theoretical models to explain 
this behavior; however, we show that their observations 
are fully consistent with the simpler YRZ model. 

In Section 2 we summarize the basic elements of the 
model that arc needed to calculate the AC optical con- 
ductivity. In Section 3 we present formulas for AC con- 
ductivity at any temperature in the pseudogap state, and 
we present numerical results for a number of dopings. Wc 
focus on the approach to the Mott insulating state, charac- 
terizing the elements that control the size of the remaining 
Drude response: the reduction in number of charge carri- 
ers, the reduction of the coherent response due to correla- 
tions, and the reduction in carrier mobility. In Section 4 we 
consider optical properties in the superconducting state. 
We examine the effect of pseudogap formation on the ther- 
mal conductivity and low-frequency optical conductivity, 
showing modifications to the Wiedemann- Franz law at sig- 
nificant temperatures but no effect in the approach to the 
universal zero-temperature limits. We then consider the 
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microwave conductivity, comparing numerical results to 
experimental data on ortho-I and -II YBCO; we find that 
the YRZ model provides qualitative agreement with the 
data, while Fermi liquid theory does not. In Section 5 we 
summarize our findings. 



2 Formalism 

In the phenomenological resonating valence bond spin liq- 
uid model of Yang, Rice, and Zhang, the coherent part of 
the charge carrier's Green's function is written as 



9t{x)Wg 



where the two energy branches are given by 



with 



+ Z\2 



The two branches have weights 
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(5) 



In these formulas, k is momentum and ui is energy; Apg 
is the pseudogap; 

^k = — 2t{x){cos kxtt + cos kyu) — 4:t'{x) cos kxa cos kyU 
— 2t"{x) {cos 2kxa + cos 2kya) — fip (6) 

is the band-theory dispersion curve for a tight-binding 
Hamiltonian including contributions up to third-nearest- 
neighbor hopping, with a chemical potential /Uj,; and = 
— 2t(a:)(cos kxa+cos kyO) is the first-nearest-neighbor con- 
tribution alone. The hopping parameters are t{x) = gt{x)to 
+ l9s{x)Jx, t'{x) = gt{x)t'Q, and t"{x) = gt{x)tQ, where 



2x 4 



(7) 



are Gutzwiller factors accounting for the loss of coher- 
ence due to correlations and narrowing of bands. The pa- 
rameters J and X take values J = |fo and x = 0.338, 
and in accordance with YRZ, we use I'q — —0.3to, and 
tQ = 0.2^0, chosen to duplicate the energy dispersion of 
Ca2Cu02Cl2. Both the pseudogap and the superconduct- 
ing gap are taken to have d-wave symmetry. 



Apg{x){cos k^a — cos ky a), (8) 
A° ^{x){cos kxtt — cos kyo), (9) 



with doping-dependent magnitudes 

A°g{x) = 0.6to(l - x/0.2), (10) 
A%{x) = O.Uto[l - 82.6{x - 0.2)2]. ^^^^ 

These dependencies place both optimal doping and the 

quantum critical point for pseudogap formation at a; = 
0.2. One could easily alter this placement, but our aim is 
to analyze qualitative features rather than fit the model 
to experiment. Superconductivity has been included at 
the level of BCS theory modified to account for rf-wave 
gap symmetry, and we give the gap a temperature depen- 
dence to accord with that. We take the pseudogap to be 
independent of temperature. 

One can use the Green's function (1) to find the asso- 
ciated regular and anomalous spectral functions 



A{k,co) = J2 9tW^[{u"f6iuj - E^) 



a=± 



+ {v''y6{co + E^) 
B{k,co)^J29tW;:^[S{ 
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E^) - S{co + E^)], (13) 
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Here we see that there arc four energy branches, iEg, 
where Eg = ^ {E^) + A'j^. To account for impurity scat- 
tering, in all our calculations we broaden the Dirac delta 
functions in the spectral functions into a Lorentzian form 
with half-width r/2. 

Many of our results in this paper can be understood 
from the manner in which the pseudogap in this model 
causes the large Fermi surface of Fermi liquid theory to 
reconstruct into a hole Luttingcr pocket centered on the 
nodal direction and an electron pocket centered on the 
antinodal direction in the Brillouin zone. Each of these 
pockets has a strongly weighted side, which follows fairly 
closely the underlying Fermi liquid Fermi surface, and a 
weakly weighted side, which is more closely related to 
the antiferromagnetic Brillouin zone boundary defined by 
^£ = 0. As the doping is decreased, these pockets shrink; 
for sufficiently large values of the pseudogap (specifically, 
for values of x < 0.16), the electron pockets are gapped 
out entirely. This situation is shown for the first quadrant 
of the Brillouin zone in the insets of the left-hand col- 
umn of Fig. 1: at X = 0.2, there is no pseudogap, and the 
Fermi surface is the large open contour of a Fermi liquid; 
at smaller values, the Fermi surface is reconstructed into 
ever-smaller pockets. But though these pockets shrink, 
they are present for all nonzero values of x. And on parts of 
the contours around these pockets, charge carriers main- 
tain high mobility and low-energy excitations can easily 
occur. Hence, the hole contours always provide a source of 
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mobile charge carriers and maintain coherence and metal- 
licity in the system as the Mott insulating state is ap- 
proached. As we shall see in the next section, this leads to 
a persistent Drudc response in the AC optical conductiv- 
ity. In addition, the reconstruction of the Fermi surface 
leaves the area around the node unaffected, preserving 
low-temperature, low-energy physics that accesses states 
only near the node. In Sec. 4 we show that this leads to 
traditional zero-temperature universal limits persisting in 
the pseudogap state. 

For om- purposes, the Green's function in Eq. (1) can 
be viewed as phenomenogical, with its primary use be- 
ing a qualitative understanding of the anomalous proper- 
ties associated with pseudogap phenomena in underdoped 
cuprates. However, its form is based on a microscopic the- 
ory of arrays of two-legged Hubbard ladders in a doped 
spin liquid including long-range interladder hopping [1]. 
And although it was initially considered only for weak 
coupling, comparison with the numerical results of Troyer 
et al. [19] for strong coupling t-J ladders suggests that it 
is characteristic of all doped spin liquids. Therefore, the 
physics encapsulated in Eq. (1) should apply to any ma- 
terials that can be described by the microscopic Hubbard 
model near half-filling. 



3 AC conductivity and optical effective mass 

In this section we examine the optical conductivity as the 
insulating state is approached. Since that approach lies 
outside the superconducting dome, we focus on the non- 
superconducting state, even when analyzing dopings and 
temperatures that would fall within the dome in exper- 
iment. Now, the real part of the AC conductivity a = 
f 1 -|- i(T2 at frequency u and temperature T is given by 



(7i(w,r) 



27re2 



duj' [/(u;' + a;)-/(a;')] 



A{k,uj')A{k,uj' +uj) 



+ B{k,uj')B{k,oj' + uj) , 



(16) 



where e is the magnitude of the electron's charge, Vk^ is 
the X component of the velocity in a state of momentum 
k, and f{oj) is the Fermi function l/[l-|-exp(a;/A:BT)]. The 
sum over k extends over the first Cu02 Brillouin zone. In 
this expression, we have neglected vertex corrections: the 
only effect of these is to change the overall magnitude of 
the conductivity, which is not essential to this work. 

In the non-supcrconducting state, the imaginary part 
of (7 can be obtained from the real part via a Kramers- 
Kronig transformation. In the superconducting state, the 
formation of the superfluid transfers part of the conduc- 
tivity into a Dirac delta function at w = 0. Including the 
Kramers-Kronig transformation of this contribution leads 
to 

a2{uj,T) = KK[a^]{uj,T) + ^^, (17) 



where -fsr-R'[tTi] is the Kramers-Kronig transformation of 
Eq. (16) and Wgc is the optical weight transferred into 
the delta function at zero frequency. Wgc can be obtained 
from the sum rule 



2 



0+ 



af)duj, 



(18) 



where ct^ and erf are the optical conductivities in the nor- 
mal and superconducting states, respectively. This stan- 
dard rule from BCS theory has been shown to hold in the 
YRZ model as well [5]. 

As we can see in Eq. (16), the conductivity depends 
on the product of two spectral functions, A(fc,w') and 
A{k, Lo'+uj). In the superconducting state, there is also the 
similar product of Gorkov anomalous functions B{k,uj') 
and B{k,uj' + lj). These products describe the transition 
from energy w' to lu' plus the photon energy uj. Each spec- 
tral function describes two branches, a = + or — , with re- 
spective weights and Wj^ . The product of two spectral 
functions hence contains four terms, which we will denote 

by ++, , H — , and — +. The first two of these, -|--|- 

and , describe intraband transitions within the + ov — 

band; the latter two describe interband transitions. 

We can easily determine that the intraband transitions 
create a Drude response centered at w = 0, while the 
interband transitions create a peak centered at some finite 
frequency. Evaluating Eq. (16) in the clean limit (r" — ^ 0), 
we find (Ti = an + criB, where 



ao = -2Tre^ 9t 



criB 



k a=± 
,.2 TJA+T 



S{uj), (19) 



X [5{uj -E+- Es) + Sioj + E+ + Es)] 
Es - Eg 

X [5{uj -E+ + Eg) + 6{u; + E+ - E^)] I. (20) 



The term au, peaked at w = 0, is the Drudc response 
coming from intraband transitions, while the term ct/b, 
with contributions at finite energies ±Eg ± Eg , is the 
interband contribution arising from transitions between 
±Ej and ±Eg . In the non-superconducting case, we can 
reproduce the exact numerical results with impurities (i.e. 
those produced with Lorentzians of half-width r/2 in the 
spectral functions) by broadening the delta functions in 
the clean-limit formulas into Lorentzians of half-width F. 
(One should note that in the superconducting case at low 
temperature, this procedure fails, since there is a discon- 
tinuity at r = 0. The difference between the supercon- 
ducting and normal states can be understood from the 

—d{Eg); this delta function, 



fact that at T = 0, 



4 



Adam Pound et al.: Optical properties of the pseudogap state in underdoped cuprates 




800 



Fig. 1. (Color online) The optical conductivity at various dop- 
ings in the non-superconducting pseudogap state at tempera- 
ture T = O.Olto/fcs and with impurity scattering F — O.lto- 
Left column: the real part cti (u) of the conductivity as a func- 
tion of frequency. The insets show the Fermi surface in the 
upper-right quadrant of the first Brillouin zone; the dashed 
green line indicates the boundary of the antiferromagnetic Bril- 
louin zone. Middle column: the corresponding imaginary part 
(72 (i^). Right column: the partial optical sum N{ljj), normal- 
ized to one-eighth the square of the plasma frequency. In each 
plot, the red shaded region indicates the contribution from in- 
traband scattering; the blue, from interband. 
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Fig. 2. (Color online) Experimental results for optical conduc- 
tivity in LSCO (left panels) and YBCO (right panels) with a 
sequence of dopings. Red shaded regions indicate Drude con- 
tributions, while blue regions indicate contributions from mid- 
infrared Lorentzian oscillators. The top two sets of panels are 
for nonsuperconducting crystals at lOK. The other panels are 
for superconducting crystals at T ~ Tc (or T* , the tempera- 
ture at which the pseudogap opens, for y = 6.75). This figure 
is taken from Lee et al. [18] with permission from the authors. 



which arises only in the limit _r — 0, has support on a 
surface if Age = but only at the nodal point if Age ^ 0.) 

In the left-hand column of Fig. 1, we display numer- 
ical results for cti as a function of frequency in the non- 
superconducting state. (Because we have restricted our- 
selves to a Cu02 plane, the conductivity takes on units 
Ig, where d is the distance between Cu02 planes.) We fo- 
cus on the behavior as doping is decreased. At a; = 0.20 
there is no pseudogap, and the Fermi surface (shown in 
the inset) is the familiar surface of Fermi liquid theory. 
In this case there is only one band, so only intraband 
transitions contribute to the conductivity. And the con- 
ductivity profile, shown in black, closely follows the well- 
known Drude form, with only slight modifications due to 
our model's more complicated band structure, which is re- 
stricted to the Cu02 Brillouin zone rather than the usual 
continuum approximation with infinite bands. At dopings 
smaller than x — 0.20, there is a finite pseudogap, and the 
Fermi surface is reconstructed into contours defining Lut- 
tinger hole pockets (again shown in the insets). In these 
cases the two bands a — zL appear, leading to both an 
intraband contribution (indicated by the red shaded re- 
gions) and an interband one (indicated by the blue shaded 
regions). 

Analogous sequences of results for a range of dopings 
were obtained experimentally by W.J. Padilla et al. [17] 



and Y.S. Lee et al. [18]. We show the experimental data 
for LSCO and YBCO in Fig. 2, reproduced from Ref. [18]. 
In that figure, the blue shaded regions indicate all sources 
of mid-infrared absorption, including incoherent scatter- 
ing and possibly charge- or spin-density waves, not only 
the interband contribution of our model. However, the red 
shaded regions unambiguously show the Drude conductiv- 
ity, the magnitude of which is a measure of the remain- 
ing coherence in the system as the doping is decreased 
toward the Mott insulating state at x = 0. Our results 
are in agreement with the data, showing the Drude re- 
sponse greatly diminishing with decreasing x (note the 
changing scale on the y-axis in both figures), but remain- 
ing finite so long as x is finite. One can easily understand 
this within our model: As x is decreased, the overall scale 
of the conductivity is greatly decreased, primarily due to 
the Gutzwiller factor gf appearing in Eqs. (19) and (20) 
but also due to the decreasing size of the Luttinger hole 
pocket, which reduces the size of the Fermi surface and the 
number of available charge carriers. But the hole pocket 
remains of finite size for a: > 0, which allows the coherent 
response to persist. 

In the middle column of Fig. 1, we display the corre- 
sponding imaginary part of the conductivity. In the pseu- 
dogap regime there are two distinct contributions to this 
quantity as well. The intraband contribution is as ex- 
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pected for a Dniclc response, but the interband piece is 
quite different in form, becoming negative at small uj. 

In the right-hand column, we display the partial opti- 
cal sum 



iV( 



Lj) = / ai{ij')dui'. (21) 
Jo 

We have normalized this quantity to the total spectral 
weight ujp/8 = Gi{Lo')du]' , where u)p is the plasma fre- 
quency. For the Fermi liquid case (x — 0.2), a single en- 
ergy scale is present, while in the pseudogap state a second 
scale arises, indicated by the small humps in the curves. 
From the blue and red shaded regions, we see that the 
interband transitions contribute a larger fraction of the 
optical weight as the doping is decreased; eventually, as 
shown in the upper-right-hand frame, the interband con- 
tribution dominates over the intraband one. Nevertheless, 
in the YRZ model, the Drude response persists for all pos- 
itive values of doping. 

We can determine more about the approach to the in- 
sulating state by examining two additional quantities: the 
effective optical mass TO*p of the charge carriers, and the 
effective number n^a of them. The effective optical mass 
is defined by a generalized Drude form, which has become 
fairly standard in disciissions of the optical conductivity 
of highly correlated systems. We write the conductivity in 
terms of m*p and an optical scattering rate 1/r: 



Air \/t{uj) — iu)m*p{u))/m 



(22) 



where m is the bare band mass. From this formula, we see 
that m*p is given by 



"op 



[CO) 



1 



0-2 



TO 



(23) 



Here TO*p is a function of frequency, but we will be in- 
terested only in the limit Y\m^^Qm^^{(jj) /m] henceforth, 
m*p/m will denote this limiting value. 

Our numerical results for function of doping 

are shown in the lower frame of Fig. 3. The mass depends 
on frequency, temperature, impurity scattering, and the 
state of the system (i.e. whether or not a pseudogap is 
present, and whether or not the state is superconducting). 
When the system is non-superconducting, the effective 
mass becomes equal to the bare mass at a; = 0.2, where 
the pseudogap vanishes. When a pseudogap is present at 
lower dopings, the effective mass increases, but it changes 
by only a factor of 1-4 across the phase diagram. Fur- 
thermore, this increase is suppressed by impurity scatter- 
ing: when r = O.lto and F = 0.2to, the effective mass is 
roughly constant, and clearly will not diverge in the ap- 
proach to a; = 0. (Temperature causes an analogous sup- 
pression, but only at the experimentally irrelevant values 
T>0.1to/kB.) 

Although our primary interest in TO*p is its behavior 
in the non-superconducting state, we also present results 
when both gaps are present (dotted curve in Fig. 3). In 



this case the mass diverges at the bottom of the supercon- 
ducting dome, at a; « 0.9 in our phase diagram. Also, the 
effective mass does not equal the bare mass at a; = 0.2. 
However, when the system is superconducting, the effec- 
tive mass defined by (23) no longer has a straightforward 
interpretation as an electron mass. Instead, it is related 
to the superfluid density. Specifically, it is the ratio of to- 
tal optical weight to the optical weight contributed by the 
superfluid: 



op 
TO 



(24) 



which follows immediately from Eq. (23), given that the 
term in Eq. (17) blows up in the limit w — >■ 0. At 
the bottom of the superconducting dome, the superfluid 
density vanishes, causing the right-hand side of Eq. (24) 

to diverge. 

We now define an effective number of charge carriers 
involved in the Drude response, beginning with the stan- 

dard formula n = . Because we are interested only in 
the carriers involved in the Drude response, we replace 
with AWd, where Wd is the Drude optical weight defined 
by 

-^=J^ <yD{^)d^. (25) 

The conductivity contains a factor (a;), where gt{x) gives 
the quasiparticle weight remaining in the electronic Green's 
function, with the rest transferred by correlations to an 
incoherent background. As we wish to count the effective 
number of carriers, we divide this factor out of our def- 
inition of rieff. Furthermore, we wish to account for the 
changes of the effective mass away from the band mass, 
so we use to.* in place of m. Our definition is thus 



rieff 



1 -n^l^Wp 
.9? 7re2 



(26) 



The dashed red curve in the upper frame of Figure 3 shows 
rieff as a function of doping. For comparison, we also show 
(in solid black) the result when the band mass m, rather 
than TO*p, is used in the definition of rieff' ■ In both cases, 
the defined number of charge carriers decreases toward 
zero in the approach to the insulating state. (Though the 
dashed curve appears to be approaching a finite value, this 
is a misleading artefact of the relatively large increase in 
effective mass in the clean limit; since the two curves differ 
by a factor of the effective mass, so long as the effective 
mass does not diverge, the dashed curve must approach 
zero if the solid curve does.) Furthermore, below x « 0.16, 
where only the hole pockets exist, the effective number 
of charge carriers is directly proportional to the doping, 
confirming that our definition of ricff is appropriate. Above 
X K, 0.16, the behavior changes as the electron pockets 
appear together with the hole pockets. 

We can now clearly understand the essential features of 
the optical response in deeply underdoped cuprates. The 
Drude conductivity decreases rapidly as a; — )■ 0, but it re- 
mains finite for all positive values of a;. In the foregoing 
discussion, we have identified three factors governing this 
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Fig. 3. (Color online) Top: the effective number Ueff of charge 
carriers as a function of doping x in the non-supercondicting 
state with no impurity scattering. For the dashed red curve, the 
renormalized optical mass is used to define riofi; for the solid 
black curve, the bare mass is used. Bottom: the optical effective 
mass m*p as a function of doping, normalized to the bare mass 
m. Note the dependence on the residual impurity scattering 
rate F (expressed here in units of to) - The dotted curve is in the 
superconducting state in the clean limit, and it diverges at the 
end of the superconducting dome, corresponding to a:: ~ 0.09 
in our model's phase diagram. All curves in both panels are for 
zero temperature. 




Fig. 4. (Color online) Conductivity and specific heat in the su- 
perconducting state, with a pseudogap (solid curves) and with- 
out (dashed curves) . Dotted red and dot-dashed blue curves are 
in the non-superconducting state. Shaded regions emphasize 
the difference between the pseudogap and Fermi liquid results. 
Top left: the DC electrical conductivity ctq = ai{uj = 0) as 
a function of temperature. Top right: the corresponding ratio 
of thermal conductivity k/T to the DC conductivity. Lower 
left: the specific heat 7 as a function of temperature. Lower 
right: the real part of the optical conductivity as a function of 
frequency. All curves are for a; = 0.12 and F — O.Olto 



behavior. First, the square of the Gutzwiller factor gtix) 
tells us that the remaining weight in the quasiparticle 
pole decreases as correlations grow and transfer spectral 
weight to the incoherent background. Second, the number 
of charge carriers drops linearly with decreasing doping. 
Third, the optical effective mass increases. However, only 
the first two factors have a large impact on the conduc- 
tivity; the increased mass has only a modest impact, since 
the carriers become heavier by only a factor of roughly 
1-4 at dopings as low as 1%. 



4 DC conductivity, Wiedemann-Franz law, 
and microwave conductivity 

In the preceding section we focused on the non-super- 
conducting state. When a superconducting gap is included, 
all points on the Luttinger pocket contours become gapped 
except those in the nodal direction. But from Eq. (5) we 
see that the weight associated with the backside of 
the pocket (i.e. the side facing the antiferromagnetic Bril- 
louin zone boundary) is zero in the nodal direction. Hence, 
there is effectively only one point that remains ungapped: 
the point in the nodal direction on the side of the pocket 
facing toward the center of the Brillouin zone. This is re- 
ferred to as the Dirac point. At very low temperatures and 
photon frequencies, this is the only active point. But this 
point and those in its very near vicinity are largely unaf- 
fected by pseudogap formation, since the pseudogap van- 



ishes in the nodal direction. Therefore, even though the 
pseudogap radically reconstructs the Fermi surface into 
small pockets, at sufficiently low temperatures it has neg- 
ligible effect on low-frequency transport properties such as 
the DC electrical conductivity, thermal conductivity, and 
microwave conductivity. In this section we show the valid- 
ity of these statements and the effect of the pseudogap as 
the temperature is increased. 

We obtain the DC conductivity by letting w — >■ in 
(16), yielding 



-< 



duj' 



, df{uj') 
duj' 



[A{k,iu'f + B{k,uj'f] , 



(27) 



where do = limt^_>.o (Ti(w). The thermal conductivity is 
given by the similar formula 



T 



J^2 



duj' 



(28) 



The top-left frame of Fig. 4 shows cto as a function of T. 
We consider two cases: the results of the YRZ model with 
a pseudogap (solid curve) and the results with Apg set to 
zero (dashed curve). As we expect, at low temperatures 
the two curves are indistinguishable, since the pseudogap 
has not changed the environment about the Dirac point. 
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As the temperature increases, however, the DC conduc- 
tivity in the superconducting state is suppressed below its 
Fermi hquid value, since the pseudogap gaps out states 
that would otherwise be available. The shaded region in 
the figure emphasizes this effect. Since the thermal con- 
ductivity exhibits the same behavior, we do not explicitly 
show it. Instead, in the upper-right frame, we show the 
Wiedemann- Franz ratio . Again the shaded region em- 
phasizes the diff'erence between results with and without 
a pseudogap, and again, at low temperatures the pseu- 
dogap is seen to have no effect. Analogous results were 
found by Borne et al. [9] in their study of specific heat 
7 = cy/T, and for comparison we include those results in 
the lower-left frame. 

As mentioned previously, we expect the same behav- 
ior in the case of the microwave conductivity: given a low 
temperature, if the photon frequency is sufficiently small, 
the pseudogap should have a negligible effect. Our expec- 
tations are confirmed by a plot of cti as a function of 
frequency at zero temperature, shown in the lower-right 
frame of Fig. 4. As we expect, the curve in the Fermi 
liquid case (dashed) merges with the curve in the pseu- 
dogap case (solid) at low frequencies. For the case shown, 
the pseudogap becomes important only in the frequency 
range where interband transitions predominate; for the 
Fermi liquid case, in which no such transitions occur, the 
conductivity behaves simply as a Drude tail in this region. 

Figure 5 displays results for the very-low-temperature 
behavior of electrical and thermal conductivity. In a pre- 
vious paper, Carbotte [20] derived an analytic expres- 
sion for the universal (i.e. impurity-independent) zero- 
temperature limits of both thermal and electrical conduc- 
tivity in a d-wave superconductor with a pseudogap. These 
limits are important signatures of nodes in the supercon- 
ducting order parameter [21,22]. Including the leading- 
order temperature correction (which shows the first de- 
pendence on impurity content), the limiting behavior is 
given by 

T^^*T^(,^+15^ 7^ + ---j' (30) 

where vp = and va = | ^dk \vd respectively 

the Fermi and gap velocities at the Dirac point Vdi a-nd 
we have set h = d = ks = ^, except in the overall pref- 
actor of fc^. Except for the important factor of gj, these 
formulas are unaltered from standard results for d-wave 
superconductors [23-25]. They are derived by performing 
a Sommerfeld expansion and expanding around Vd- In 
that derivation, one must assume the residual impurity 
scattering rate F is much smaller than all energy scales 
other than T; for small scattering rates, one can show 
that Eqs. (29) and (30) are accurate to within an error of 
order . We see in Fig. 5 that the analytic expressions 
yield excellent agreement with numerical evaluations of 
the exact formulas. As we have shown above, because the 
pseudogap shares the c?-wave symmetry of the supercon- 
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Fig. 5. (Color online) DC electrical conductivity ctq (open blue 
triangles) and thermal conductivity k/T (open black circles) as 
a function of temperature (in units of residual scattering rate 
r) in the superconducting state. The dashed and dotted red 
lines are analytic forms valid in the limit T/ F — >■ 0. All curves 
are for a; = 0.12 and F = 0.005to- 

ducting gap, these results are valid regardless of whether 
or not the pseudogap is present. Note that these results 
hold only for finite F; as mentioned in the preceding sec- 
tion, in the superconducting state there is a discontinuity 
at (7^ = 0,T = 0), and if F vanishes identically, the con- 
ductivity at r = also vanishes identically, rather than 
reaching the universal value given above. 

All derivations of these limits involve only the coher- 
ent part of the electron propagator. However, Schachinger 
and Carbotte [26] considered the case of coupling to a bo- 
son spectrum within Eliashberg theory. Such a calculation 
fully accounts for the inelastic scattering that provides the 
boson-assisted incoherent contribution to the conductiv- 
ity. They showed analytically that in this case Eqs. (29) 
and (30) still hold, but with the factor of gf replaced by 
1/(1 -I- A), where A is the mass-enhancement factor asso- 
ciated with the electron-boson interaction. Like g^, the 
factor 1/(1 -I- A) measures the reduction of the coherent 
part of the conductivity due to interactions. 

We now move to our final results, shown in Fig. 6. 
These are again plots of the DC conductivity as a func- 
tion of temperature, but unlike in the preceding plots, 
which used a constant impurity scattering rate, here we 
have used a temperature-dependent rate in order to ac- 
count for inelastic scattering. This allows us to meaning- 
fully compare with experimental data on microwave con- 
ductivity: it is now widely accepted that the large peak, 
usually at a temperature much below Tc, in the microwave 
conductivity [28-31] and thermal conductivity [32,33] of 
cuprates is not a coherence peak associated with super- 
conductivity, but rather the effect of a rapid decrease in 
inelastic scattering with decreasing temperature below Tc- 
One prominent model that has emerged from fits to data 
is a power law [31,34], and here we have used 

r(T) = O.OOKo + Q.mQ{T/Tcf (31) 

for all values of doping. In the upper frame of Fig. 6, we 
show results for Fermi liquid theory, with no pseudogap 
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Fig. 6. (Color online) DC conductivity ctq in the supercon- 
ducting state as a function of temperature for four values 
of doping X, with a temperature dependent scattering rate 
r = O.OOlto + Q.l&to{T/Tcf. Top: results in the Fermi liq- 
uid model. Bottom: results in the YRZ model. Results in the 
Fermi liquid model are related to those in the YRZ model by 
setting the pseudogap to zero and dropping all Gutzwiller co- 
herence factors (except those that define the band structure). 
The inset in the bottom frame shows experimental microwave 
conductivity results for ortho-I and -II YBCO, taken from R. 
Harris et al. [27]. For ortho-I, two samples were considered, 
with critical temperatures Tc = 89K and 91K. For ortho-II, 
Tc = 56K. All experimental data is for a-axis conductivity at 
l.lGHz. 



present and no overall factor of in the conductivity. 
(However, the Gutzwiller factors defining the band struc- 
ture are left intact, and we have scaled the vertical axis 
by the Gutzwiller factor at optimal doping to allow easy 
comparison with the YRZ model.) In the lower frame, we 
show results for the YRZ model. In both cases, we show 
a range of dopings, and we see that though the shapes of 
the curves are very similar in both models, the Gutzwiller 
factor in the YRZ model causes the magnitude of the con- 
ductivity to decrease much more rapidly with decreasing 
doping. (We have confirmed numerically that the pseudo- 
gap has only a negligible effect in comparison.) 

The inset in the lower frame shows experimental data 
for microwave (1.1 GHz) conductivity, taken from Harris et 



al. [27]. The upper data set was obtained from two samples 
of optimally doped ortho-I YBa2Cu30e.g93 with T^, values 
of 89K and 9 IK; the lower set, from a sample of ortho-II 
YBa2Cu306.5 with Tc = 56K. For ortho-II, every second 
chain is empty. Given our choice of a; = 0.2 as optimal 
doping, ortho-II would correspond to a doping between 
X = 0.10 and x — 0.13. Comparing these data sets with 
the corresponding numerical curves, we see that the YRZ 
model, with its large decreases in magnitude with decreas- 
ing doping, yields far better qualitative agreement with 
experiment than does the Fermi liquid model. In order to 
match the experimental data, we could alter the param- 
eters in the scattering rate, and in particular, we could 
allow them to differ between ortho-I and ortho-II. (We 
could also generate numerical results at l.lGHz, rather 
than DC data, which would slightly lower our curves.) 
However, here we are interested only in qualitative behav- 
ior, and our results clearly show that the data is more eas- 
ily understood within the YRZ model than ordinary Fermi 
liquid theory. Furthermore, we see that the Gutzwiller fac- 
tor is the essential element in explaining the data. 



5 Summary and conclusions 

We have found that the approach to the insulating Mott 
state captured in the resonating valence bond spin liquid 
model developed by Yang, Rice, and Zhang is in quali- 
tative agreement with the experimental observation that 
a coherent Drude response persists in the deeply under- 
doped regime of the cuprates down to the 1% doping level. 
In this regime, small hole Luttinger pockets remain, and 
on the boundary contours of these pockets zero energy 
excitations exist, leading to a Drude term in the conduc- 
tivity. This term decreases as the size of the pockets (and 
therefore the number of charge carriers) shrinks with de- 
creasing doping. It also decreases due to the presence of 
Gutzwiller factors accounting for increased correlations, 
which cause a transfer of spectral weight to an incoher- 
ent background. An increase in optical effective mass also 
causes a decrease in metallicity as the insulating state is 
approached, but in agreement with experiment, this is not 
a dominant effect: in the regime of interest, the effective 
mass is only of order twice the band mass. In this sense, 
the charge carriers do not become heavy. 

We have also investigated the effect of pseudogap for- 
mation on the temperature dependence of the DC and 
thermal conductivity in the superconducting state. If the 
Gutzwiller factor is not considered, the pseudogap leaves 
unaltered the low-temperature behavior of the conductiv- 
ities, even though it causes a radical reconstruction of the 
Fermi surface, from the large open contour of Fermi liq- 
uid theory to small Luttinger pockets. This observation 
is explained by the fact that the zero-temperature con- 
ductivities depend only on the region of reciprocal space 
around the tip of the Dirac cone, which is unaffected by 
the finite pseudogap because the pseudogap has d-wave 
symmetry, vanishing along the nodal direction where the 
Dirac point lies. Due to this, the universal limit associated 
with d-wave superconductivity survives in the pseudogap 
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regime. However, the universal limit is now modified by a 
Gutzwiller factor (a;). The same holds true for the ap- 
proach to zero temperature: analytic formulas derived for 
the small-T limit are unaltered except for the appearance 
of a Gutzwiller factor. 

Finally, we have considered microwave conductivity. 
Plots of conductivity as a function of temperature show 
that Fermi liquid theory and the YRZ model yield curves 
of very similar shape but greatly different magnitude. Com- 
paring these numerical results to existing data for ortho-I 
and -II YBCO, we found that only the YRZ model is in 
qualitative agreement with experiment. The lower conduc- 
tivity of ortho-II as compared to ortho-I is seen to result 
primarily from the increased correlations described by the 
Gutzwiller factor. 

We wish to thank D. N. Basov and Y. S. Lee for allowing us 

use of a figure from Ref. [18] . This work was supported in part 
by the Natural Sciences and Engineering Research Council of 
Canada and the Canadian Institute for Advanced Research. 
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